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A NOTE ON THE PLASMA MOMENT 
EQUATIONS IN A DIPOLE FIELD 

J. M. Grebowsky 
ABSTRACT 

A coordinate system is discussed which is defined with respect 
to a magnetic dipole field by unit vectors in the direction of the field, 
in the azimuthal direction and in a direction along the dipole equipo- 
tential lines. The velocity terms of the continuity and momentum 
equations pertinent to plasma motion are expressed in terms of vari- 
ations in these directions. 
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A NOTE ON THE PLASMA MOMENT 
EQUATIONS IN A DIPOLE FIELD 

1. Introduction 

The magnetic dipole field plays an important role in determining the prop- 
erties of the plasma distribution and the nature of the waves propagating within 
the plasmasphere since the magnetic field in the region is approximately dipolar. 
Therefore coordinates should be used which reflect the geometry of the field 
keeping in mind that an investigation of the plasma moment equations is simpli- 
fied by expressing vector quantities in terms of components along the directions 
in which the principal processes take place. 

Gothard (1967) has summarized the different descriptions of the dipole field 
and has investigated the mathematical form of the differential distances along 
and across the field lines. This paper will develop the mathematical forms of 
the continuity and momentum equations in a curvilinear coordinate system de- 
fined at every point by unit vectors in the direction of the field (i„ ), in the 
direction perpendicular to the field direction but lying in the plane of the field 
line ( i| ), and in the azimuthal direction about the dipole axis ( i^). The orienta- 
tions of these vectors are depicted in Figure 1. 

The distance along a field line is the coordinate with respect to which the 
thermal plasma moment equations are frequently expressed in studies of 
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ionosphere-protonosphere coupling (see for example Angerami and Thomas, 
1964, and Tamao, 1966). The azimuthal direction will correspond to the direc- 
tion of the plasma velocity for corotation about the earth's dipole axis. The 
direction ij_ , which is identical to the radial direction in the equatorial plane, 
will be useful in studies of diffusion across field lines. Although many of the 
terms to be investigated have previously been used in ionospheric and magneto- 
spheric computations, their complete expression in terms of these coordinates 
has not been previously developed. 
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2 . The Coordinate System 

In terms of spherical coordinates (r, 0 , 9 ) where r is the geocentric dis- 
tance, 8 is the colatitude and 9 is the azimuthal angle, the dipole magnetic field 
is: 


„ 2 a cos 6 « a sin 8 ^ 

B = r + 8 

r 3 r 3 

or 

B = grad cos 8 ^ 
where a is the dipole moment. 



The unit vectors along and perpendicular to the field direction (Gothard, 
1967) are 
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A * 

2 cos 6 r + sin 9 9 
/l + 3 cos 2 6 
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r - 2 cos 0 <9 + sin 0 r 


l i = 


(1+3 cos 2 0) 1/2 


i<p = 9- 


These unit vectors form an orthogonal system and are related to one another by 


l l - x n x x 9‘ 


The relative orientation between the direction ij_ and the direction of the 
dipole gradient iy is given as a function of colatitude by 


iy * ij. 


i, = 


- sin 0(1+ cos 0 ) 


t(l + 3 cos 2 0) 2 + cos 2 6 sin 2 0] 1/2 [4-3 sin 2 0] 1/2 


The length coordinates corresponding to these unit vectors are easily de- 
termined. Measured from the center of the dipole axis, the distances along 
these coordinate directions are in terms of the dipole latitude \ and the geo- 
centric distance r : 


2 v' 3 cos 2 k 


'In 


(■ 


/3 + 2 


\ /3 sin \ + /l + 3 sin 2 \ 

l = I 1 im I* 

2 sin 1/2 K e-o J e V sin \ / 


- /3 sin k S' + 3 sin 2 k + 2 /3 


1/2 


dX. 


and along the azimuthal direction 


= r <p sin 6 . 
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The distance S ( , is the distance along a field line (i.e., along the curve r = 
constant ° cos 2 \ ) whereas Sj^ is the distance along an equipotential line of the 
dipole (r 2 = constant ° sin \). The reference level for is taken as 9 = 0. 

The expression for contains an integral which has not been evaluated 
analytically. However, near the equator where sin K ~ \ the integral can be 
evaluated approximately to yield 



It should be noted that near the poles (A. = ± n / 2 ) the differential distance 
dS^ approximates to r d 9 just as near the equator dS, ( ~ r d9« However, these 
coordinates should be applied with caution near the poles since is a zero vector 

at the pole and the divergence of this unit vector has a singularity at the pole. 

In a dipole meridian plane, the contours of as defined are mirror reflections 

A 

about the dipole axis and is di rected in the polar regions towards the pole. 
Hence thevectorsinthe perpendicular direction converge at the pole forming 
the singularity. 

3 . Continuity Equation 
The continuity equation has the general form 

l£ + div J = S 

O t 
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where p is the density, J is the flux density and S represents source and sink 
terms. Only the divergence term is to be considered since the other terms do 
not explicitly involve particular spatial directions. 

Expressed in terms of the three previously discussed coordinates, the 
divergence term is 


div J = div (J,, i,| + Jj ij_+ J ? i^) 


3f 


or using the relationship i~ • grad f - : 

OS 


... T 3J " 3J i , .. • . .. * 

d,vJ = 3S^ dlVl " +J l dlv ‘i + 3sr 


where use has been made of div - 0. In terms of spherical coordinates the 


differentiation operators are 


9 _ 2 cos 9 3_ sin 9 B_ 

/l + 3 cos 2 9 r /l + 3 cos 2 6 

d sin 9 3_ 2 cos 9 3_ 

/l + 3 cos 2 9 r /l + 3 cos 2 9 ^ 
and 

a _ i a 

BS^ r sin 6 3 <p 

The divergence quantities in terms of spherical coordinates are: 
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div i t | = 


3 cos 8 (3 + 5 cos 2 8 ) 


(1 + 3 cos 2 8 ) 3/ 2 


div ij^ - 


2 sin 8 


Jl sin 2 8 - cos 2 6-3 cos 4 9 \ 
r /l + 3 cos 2 6 \ (1+3 cos 2 6 ) sin 2 6 j 


The first quantity can also be written in the familiar form: 


so that 


div i 


H 


I 3B 

R 3 S, | 


aj„/B 

divO,, '„) 


It is interesting to speculate as to what function F is required in order to 
have a similar expression for the divergence of the perpendicular vector. That 
is 


div i^ 


I IL 

F BS i 


Such a functional form would be an aid in exploring plasma phenomena in which 
the perpendicular flux is important. In terms of spherical coordinates the dif- 
ferential equation for F can be written 

3F_2c<>s_S 3F = _2 f 1+ 1 cot »4,. 

Br r sin 0*9 r | i + 3 C os 2 6 J 

Solving this equation by separation of variables yields 
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F = K r m Q + 3 cos2 2 
sin 6 (cos fc?) 1/2 m + 2 

where K is an arbitrary constant and m is determined by the geometry. 

Interpreting F as the inverse of the cross sectional area of a tube defined 
by equipotential lines (i.e., lines parallel to ), in analogy to the concept of a 
magnetic flux tube, the constant m must be -4 in order to have a non-vanishing 
area at the equator. This value can also be obtained by considering the geometry 
of an equipotential tube of infinitesimal cross section at the equatorial plane — a 
procedure which illuminates some properties of this coordinate system. 

In a meridian plane the equipotential lines are of the general form as shown 
in Figure 2. In the equatorial plane the equipotential lines are radial (see 
Figure lb). Therefore the differential cross sectional area of a potential tube 
at the equator is 


(r d 9 ) (r d | X. | ) 

where | A. | is the angle between the equator and the wall of the flux tube in a 
meridian plane. But along an orthogonal line (i.e., an equipotential) 

r = M sin 1/2 | \ | 
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wftere M is a constant. Therefore near the equator 



The differential area is then 

r 4 , 

— dtp 

M 2 

and hence the value of m is -4. For this value the function F is 


KB 

r sin 9 


and the divergence of J can be put into the particularly simple form 


j- T d ^ /B ' ) B 

div J = B — — + 


r sin 6 


B 


BS„ 


r sin 6 


BS, 


^Jcp 


Near the equator in the approximation sin A. % A. and cos \ * 1 the divergence 
becomes in spherical coordinates 


3 J|' l^Jll H 2\*h 4 
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div J - 2 X ~ + — ~ + *~ Ji + Jjj + - T~ 

dr r 3 X dr - ^ " 1 - " - ^ ^ 






r Bcp 


At the equator ( X. = 0) this expression reduces to a form which is similar to that 
of the corresponding divergence in a spherically symmetric system. That is, 
at k = 0, 
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.• /t x T N _ ^ J 1 4 T 1 B Jll 1 B Jcp 

div (J r J„ . Jqp) - 37 + 7 Ji + 7 — + 7 J^r 


whereas in a spherical coordinate system 

3 j, 


div < J '' J »' V = 37 + 7 J - + ;— *7 3^ 


1 3J 6 
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4. Momentum Equation 

The momentum equation for one component of a thermal plasma has the 
general form 

* 

d v - grad P F 
dt p p 

where P is the pressure, p the density and F denotes forces due to gravity, 
electric and magnetic fields and collisional interactions. The time derivative 
is the convective acceleration. 

In terms of the length coordinates discussed in the previous sections the 
gradient operator is 


grad = i|| _ . ^ _ + _ 


This form of the operator is needed to express the convective derivative 


dv d v - , - 

— = — + V • grad V 
dt dt 


in the desired form. 
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The local variation of the velocity with time can be written 


Bv r dv H - dv ± ? dv <p 

b 7 = x " ~ + 4 i bT + l v — ‘ 


The spatial derivative is evaluated by a straightforward expansion of the vector 
components along the i (| , and directions. The three orthogonal com- 
ponents are 


f/ - dv ll dv u dv n 2 3 cos 8 

<v ' 6 " d) ,1 " = v " ^ ^ > r /T 7 w s 


3 v, ( v. | - 2 v, C °- - % ^ sin 6 (1 + cos 2 6) 
± \ 1 sin i 9 ) 


r (1 + 3 cos 2 0 ) 3/2 


r,- dv i Bv l dv l 2 ( 1-3 cos 2 9) 

l(v • grad) Vi, = v„ + v ± — + v,, — - - v„ — 

dS|| BSjp t r s ^ n ^ (l + 3 cos 2 c?) 

, / „ cos 0\ sin 6 (1 + cos 2 6) 

+ 3 v.. v„ - 2 v. — - — 

\ sin 6 > r (1 + 3 cos 2 8) 3/2 


[(▼ - g r a d ) ▼]« = v i 


Bv„ B v„ 

t- 1 - + v. — L + v m —r 


<P~ V " 3S^ + V i3S^" >B^ 


r /l + 3 cos 2 6 


, Q (1-3 cos 2 8 ) 

3 cos 8 + V. i : — 3 i 

" sin d 


The terms quadratic in v M , v^ or v^ correspond to centripetal acceleration 
components whereas those terms which consist of the product of two different 
velocity components correspond to components of the coriolis acceleration. 
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